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Abstract. We present an improved model and theory for time-causal 
and time-recursive spatio-temporal receptive fields, obtained by a combi¬ 
nation of Gaussian receptive fields over the spatial domain and first-order 
integrators or equivalently truncated exponential filters coupled in cas¬ 
cade over the temporal domain. Compared to previous spatio-temporal 
scale-space formulations in terms of non-enhancement of local extrema or 
scale invariance, these receptive fields are based on different scale-space 
axiomatics over time by ensuring non-creation of new local extrema or 
zero-crossings with increasing temporal scale. Specifically, extensions are 
presented about parameterizing the intermediate temporal scale levels, 
analysing the resulting temporal dynamics and transferring the theory 
to a discrete implementation in terms of recursive filters over time. 


1 Introduction 

Spatio-temporal receptive fields constitute an essential concept in biological vi¬ 
sion (Hubei and Wiesel [T]; DeAngelis et al. m) and for expressing computer 
vision methods on video data (Adelson and Bergen [3]; Zelnik-Manor and Irani 
[5]; Laptev and Lindeberg [5]; Jhuang et al. [7]; Shabani et al. 0). For off-line 
processing of pre-recorded video, non-causal Gaussian or Gabor-based spatio- 
temporal receptive fields may in some cases be sufficient. When operating on 
video data in a real-time setting or when modelling biological vision computa¬ 
tionally, one does however need to take into explicit account the fact that the 
future cannot be accessed and that the underlying spatio-temporal receptive 
fields must be time-causal. For computational efficiency and for keeping down 
memory requirements, it is also desirable that the computations should be time- 
recursive, so that it is sufficient to keep a limited memory of the past that can 
be recursively updated over time. 

The subject of this article is to present an improved scale-space model for 
spatio-temporal receptive fields based on time-causal temporal scale-space ker¬ 
nels in terms of first-order integrators coupled in cascade, which can also be 
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transferred to a discrete implementation in terms of recursive filters. The model 
builds on previous work by (Fleet and Langley [9]; Lindeberg and Fagerstrom 
cni; Lindeberg |11I12] ') which is here complemented by a better design for the 
degrees of freedom in the choice of time constants for intermediate temporal 
scale levels, an analysis of the resulting temporal response dynamics and details 
for discrete implementation in a spatio-temporal visual front-end. 


2 Spatio-temporal receptive fields 

The theoretical structure that we start from is a general result from axiomatic 
derivations of a spatio-temporal scale-space based on assumptions of non-enhancement 
of local extrema and the existence of a continuous temporal scale parameter, 
which states that the spatio-temporal receptive fields should be based on spatio- 
temporal smoothing kernels of the form (see overviews in Lindeberg uma): 

T{xi,X 2, t; s, T, u, S) = g{xi - vit, X2 - V 2 t; s, E) h{t; r) (1) 


where 

— X = {xi,X 2 )"^ denotes the image coordinates, 

— t denotes time, 

— s denotes the spatial scale, 

— T denotes the temporal scale, 

— V = {vi,V 2 )'^ denotes a local image velocity, 

— E denotes a spatial covariance matrix determining the spatial shape of an 

affine Gaussian kernel g{x; s,E) = 2 TrsJdet ^ 

— g{xi — vit, X 2 — V 2 t; s, E) denotes a spatial affine Gaussian kernel that moves 
with image velocity v = (vi,V 2 ) in space-time and 

— h(t; t) is a temporal smoothing kernel over time. 

For simplicity, we shall here restrict the family of affine Gaussian kernels over the 
spatial domain to rotationally symmetric Gaussians of different size s, by setting 
the covariance matrix T" to a unit matrix. We shall mainly restrict ourselves to 
space-time separable receptive fields by setting the image velocity v to zero. 

A conceptual difference that we shall pursue is by relaxing the requirement 
of a continuous temporal scale parameter in the above axiomatic derivations by 
a discrete temporal scale parameter. We shall also replace the previous axiom 
about non-creation of new image structures with increasing scale in terms of non¬ 
enhancement of local extrema (which requires a continuous scale parameter) by 
the requirement that the temporal smoothing process, when seen as an operation 
along a one-dimensional temporal axis only, must not increase the number of 
local extrema or zero-crossings in the signal. Then, another family of time-causal 
scale-space kernels becomes permissible and uniquely determined, in terms of 
first-order integrators or truncated exponential filters coupled in cascade. 

The main topics of this paper are to handle the remaining degrees of free¬ 
dom resulting from this construction about: (i) choosing and parameterizing 
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the distribution of temporal scale levels, (ii) analysing the resulting temporal 
dynamics and (iii) describing how this model can be transferred to a discrete 
implementation while retaining discrete scale-space properties. 

3 Time-causal temporal scale-space 

When constructing a system for real-time processing of sensory data, a funda¬ 
mental constraint on the temporal smoothing kernels is that they have to be 
time-causal. The ad hoc solution of using a truncated symmetric filter of fi¬ 
nite temporal extent in combination with a temporal delay is not appropriate 
in a time-critical context. Because of computational and memory efficiency, the 
computations should furthermore be based on a compact temporal buffer that 
contains sufficient information for representing sensory information at multiple 
temporal scales and computing features therefrom. Corresponding requirements 
are necessary in computational modelling of biological perception. 

Time-causal scale-space kernels for pure temporal domain. Given the require¬ 
ment on temporal scale-space kernels by non-creation of local extrema over a 
pure temporal domain, truncated exponential kernels 

( JTp-tIu-k f > 0 

= ( 2 ) 

can be shown to constitute the only class of time-causal scale-space kernels 
over a continuous temporal domain in this sense (Lindeberg [13] ; Lindeberg and 
Fagerstrom nni). The Laplace transform of such a kernel is given by 

/.oo ^ ^ 

tlexpi.T Tk) — I hexpif] t^k) C ^ dt = - ^ (3) 

7t=-oo ^ + TkQ 

and coupling K such kernels in cascade leads to a composed filter 

l^composedif^ T) — ^k—l^^xp (i; Tk) (4) 

having a Laplace transform of the form 

pOO K 

^composedisii M) “ / (*A:—1 ^ ^ dt = I I - ^ . (5) 

Jt^-oo 

The composed filter has temporal mean and variance 

K K 

nT-K = = ( 6 ) 

k=l k=l 

In terms of physical models, repeated convolution with such kernels corresponds 
to coupling a series of first-order integrators with time constants fik in cascade 

dtL{t-, Tk) = — {L(t] Tk-i) - L{t; Tk)) (7) 

Pk 
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with L(t; 0) = f{t). These temporal smoothing kernels satisfy scale-space prop¬ 
erties in the sense that the number of local extrema or the number of zero- 
crossings in the temporal signal are guaranteed to not increase with the temporal 
scale. In this respect, these kernels have a desirable and well-founded smooth¬ 
ing property that can be used for defining multi-scale observations over time. A 
constraint on this type of temporal scale-space representation, however, is that 
the scale levels are required to be discrete and that the scale-space representation 
does hence not admit a continuous scale parameter. Computationally, however, 
the scale-space representation based on truncated exponential kernels can be 
highly efficient and admits for direct implementation in terms of hardware (or 
wetware) that emulates first-order integration over time, and where the temporal 
scale levels together also serve as a sufficient time-recursive memory of the past. 


h{t-, ij,,K = 7) ht{t\ ^j.,K = 7) 




^i,K = 7) 

A 


A 


h{t- A = 7, c = ^/2) ht{t-, K = 7, c = y/2)) htt{t-, K = 7,c = V2)) 

A 




h{t- K = 7,c = 2^/4) ht{t-, K = 7,c = 2^/^) K = 7,c = 2^/^)) 



Fig. 1. Equivalent kernels with temporal variance r = 1 corresponding to the 
composition of K truncated exponential kernels in cascade and their first- and 
second-order derivatives, (top row) Equal time constants /r. (middle row) Log¬ 
arithmic distribution of the scale levels for c = \/2. (bottom row) Logarithmic 
distribution for c = 2^/“^. 


When implementing this temporal scale-space concept, a set of intermediate 
scale levels has to be distributed between some minimum and maximum scale 
levels Train = Ti and Tmax = tk- Assuming that a total number of K scale levels 
is to be used, it is natural to distribute the temporal scale levels according to 
a geometric series, corresponding to a uniform distribution in units of effective 
temporal scale Te// = logr. Using such a logarithmic distribution of the temporal 
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hxtix,t-, s,r) -hxxt{x,t; s,r) 



h^x{x,t; s,T,v) -hxxx{x,t-, s,r,v) 


Tj ; 

n 




Fig. 2. Computational modelling of simple cells in VI as reported by DeAnge- 
lis et al. [5] using spatio-temporal receptive fields of the form T{x,t; s,t,v) = 
dx’^dti3g(x—vt\ s) h(t] t) according to equation Q with the temporal smoothing 
function h{t] t) modelled as a cascade of first-order integrators/truncated expo¬ 
nential kernels of the form Q and using a logarithmic distribution of the inter¬ 
mediate scale levels of the form ^. (left column) Separable receptive fields corre¬ 
sponding to mixed derivatives of first- or second-order derivatives over space with 
first-order derivatives over time, (right column) Inseparable velocity-adapted re¬ 
ceptive fields corresponding to second- or third-order derivatives over space. 
Parameter values: (a) hxt- <Xx = 0.6 degrees, at = 60 ms. (b) hxxt- '^x = 0.6 de¬ 
grees, Ut = 80 ms. (c) hxx'- fXx = 0.7 degrees, at = 50 ms, v = 0.007 degrees/ms. 
(d) hxxx- <Xx = 0.5 degrees, at = 80 ms, v = 0.004 degrees/ms. 


scale levels, the different levels in the temporal scale-space representation at 
increasing temporal scales will serve as a logarithmic memory of the past, with 
qualitative similarity to the mapping of the past onto a logarithmic time axis 
in the scale-time model by Koenderink El- If we have the freedom of choosing 
Train freely, a natural way of parameterizing these temporal scale levels using a 
distribution parameter c > 1 

Tfc = w {l<k<K) (8) 

which by equation (§ implies that time constants of the individual first-order 
integrators will be given by 


^l-K 


Ml — C y/ '^max 

/ - k—K—1 

l^k — y'^k '^k—1 — C 


\/c2 - ly/Tmax {“2 < k < K) 


(9) 

( 10 ) 


If the temporal signal is on the other hand given at some minimum temporal 

1 ' ' 

2(iC-l) 


i, we can instead determine c = 


in (8) such that 


scale level 

Ti = Train and add K—1 temporal scales with fj,k according to (101. Alternatively, 
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if one chooses a uniform distribution of the intermediate temporal scale levels 

k 

— yy '^max (H) 

K 

then the time constants are given by 

Figure shows graphs of such kernels that correspond to the same value 
of the composed variance, using either a uniform distribution or a logarithmic 
distribution of the intermediate scale levels. Generally, these kernels are highly 
asymmetric for small values of K, whereas they become gradually more sym¬ 
metric as K increases. The degree of continuity at the origin and the smoothness 
of transition phenomena increase with K such that coupling of iF > 2 kernels 
in cascade implies a C^^^-continuity of the temporal scale-space kernel. Specif¬ 
ically, the kernels based on a logarithmic distribution of the intermediate scale 
levels allow for faster temporal dynamics compared to the kernels based on a 
uniform distribution of the intermediate scale levels. 

Figure shows the result of modelling the spatio-temporal receptive fields 
of simple cells in VI in this way, using the general idealized model of spatio- 
temporal receptive fields in equation Q in combination with a temporal smooth¬ 
ing kernel obtained by coupling a set of first-order integrators/truncated expo¬ 
nential kernels in cascade. This result complements the general theoretical model 
for visual receptive fields in m with explicit modelling results and a theory for 
choosing and parameterizing the intermediate discrete temporal scale levels. 


4 Temporal dynamics of the time-causal kernels 


For the time-causal filters obtained by coupling truncated exponential kernels 
in cascade, there will be an inevitable temporal delay depending on the time 
constants /ifc of the individual filters. A straightforward way of estimating this 
delay is by using the additive property of mean values under convolution uik = 
Mfc according to (j^. In the special case when all the time constants are 
equal /i^ = a/t/AT, this measure is given by 

rriuni = ( 12 ) 


showing that the temporal delay increases if the temporal smoothing operation 
is divided into a larger number of smaller individual smoothing steps. 

In the special case when the intermediate temporal scale levels are instead 
distributed logarithmically according to (|^, with the individual time constants 
given by (|^ and (10), this measure for the temporal delay is given by 


'^log — 


C (c^ — (\/c^ — 1 -I- l) C -I- \/c2 — 1 c^) 


c — 1 


(13) 


with the limit value miog-umit = lim^_>oo rniog = ^ ^/t when the number 

of filters tends to infinity. 
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By comparing equations (12) and (13) including the limit value of the latter, 
we can specifically note that with increasing number of intermediate temporal 
scale levels, a logarithmic distribution of the intermediate scales implies shorter 
temporal delays than a uniform distribution of the intermediate scales. 

Table [B shows numerical values of these measures for different values of K 
and c. As can be seen, the logarithmic distribution of the intermediate scales 
allows for significantly faster temporal dynamics than a uniform distribution. 


Temporal mean values of time-causal kernels 


K 


miog (c = 

\/2) TTilog (c — . 

2^/^) miog (c = 2) 

2 

1.414 

1.414 

1.399 

1.366 

3 

1.732 

1.707 

1.636 

1.549 

4 

2.000 

1.914 

1.777 

1.641 

5 

2.236 

2.061 

1.860 

1.686 

6 

2.449 

2.164 

1.910 

1.709 

7 

2.646 

2.237 

1.940 

1.721 

8 

2.828 

2.289 

1.957 

1.732 


Table 1. Numerical values of the temporal delay in terms of the temporal mean 
m = Mfc in units of cr = y/r for time-causal kernels obtained by coupling K 

truncated exponential kernels in cascade in the cases of a uniform distribution of 
the intermediate temporal scale levels = kr/K or a logarithmic distribution 
Tfc = c^kk-K)T c > 1. 


Temporal delays from the maxima of time-causal kernels 


K t 


•max ,um 


tmax,log (c V^) tmax^Log (c 2 ) tmax^log (c 


2 

0.707 

0.707 

0.688 

0.640 

3 

1.154 

1.122 

1.027 

0.909 

4 

1.500 

1.385 

1.199 

1.014 

5 

1.789 

1.556 

1.289 

1.060 

6 

2.041 

1.669 

1.340 

1.083 

7 

2.268 

1.745 

1.370 

1.095 

8 

2.475 

1.797 

1.388 

1.100 


Table 2. Numerical values for the temporal delay of the local maximum in units 
oi a — y/r ior time-causal kernels obtained by coupling K truncated exponential 
kernels in cascade in the cases of a uniform distribution of the intermediate 
temporal scale levels = kr/K or a logarithmic distribution 
with c > 1. 


Additional temporal characteristics. Because of the asymmetric tails of the time- 
causal temporal smoothing kernels, temporal delay estimation by the mean value 
may however lead to substantial overestimates compared to e.g. the position 
of the local maximum. To provide more precise characteristics in the case of 
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a uniform distribution of the intermediate temporal scale levels, for which a 
compact closed form expression is available for the composed kernel 


hcomposed{i: K'j 


^K-1 g-t/M 


(14) 


let us differentiate this function dt (hcompoaed{t\ 
and solve for the positions of the local maximum 


e A* 

(3 r{K) 


^max.uni 




(15) 


Table shows numerical values for the position of the local maximum for both 
types of time-causal kernels. As can be seen from the data, the temporal response 
properties are significantly faster for a logarithmic distribution of the intermedi¬ 
ate scale levels compared to a uniform distribution and the difference increases 
rapidly with K. These temporal delay estimates are also significantly shorter 
than the temporal mean values, in particular for the logarithmic distribution. 

If we consider a temporal event that occurs as a step function over time {e.g. a 
new object appearing in the field of view) and if the time of this event is estimated 
from the local maximum over time in the first-order temporal derivative response, 
then the temporal variation in the response over time will be given by the shape 
of the temporal smoothing kernel. The local maximum over time will occur at a 
time delay equal to the time at which the temporal kernel has its maximum over 
time. Thus, the position of the maximum over time of the temporal smoothing 
kernel is highly relevant for quantifying the temporal response dynamics. 


5 Computational implementation 

The computational model for spatio-temporal receptive fields presented here is 
based on spatio-temporal image data that are assumed to be continuous over 
time. When implementing this model on sampled video data, the continuous 
theory must be transferred to discrete space and discrete time. 

In this section we describe how the temporal and spatio-temporal receptive 
fields can be implemented in terms of corresponding discrete scale-space kernels 
that possess scale-space properties over discrete spatio-temporal domains. 


5.1 Discrete temporal scale-space kernels based on recursive filters 

Given video data that has been sampled by some temporal frame rate r, the 
temporal scale at in the continuous model in units of seconds is first transformed 
to a variance r relative to a unit time sampling 


(16) 
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where r may typically be either 25 Hz or 50 Hz. Then, a discrete set of inter¬ 
mediate temporal scale levels is defined by ([^ or 0 with the difference 
between successive scale levels according to At^ = (with tq = 0). 

For implementing the temporal smoothing operation between two such adja¬ 
cent scale levels (with the lower level in each pair of adjacent scales referred to 
as fin and the upper level as font), we make use of a first-order recursive filter 


foutit) - foutit - 1 ) = 


1 + /tfe 


U^n{t) - foutit - 1 )) 


with generating function 


Hgeomi^) s 

1 - /ife (z - 1) 


(17) 


(18) 


which is a time-causal kernel and satisfies discrete scale-space properties of guar¬ 
anteeing that the number of local extrema or zero-crossings in the signal will 
not increase with increasing scale (Lindeberg [13) : Lindeberg and Fagerstrom 
cni). Each such filter has temporal mean value rrik = fJ'k and temporal variance 
Atj. = fff. fik, and we compute from Atj. according to 


i/r^p43rfc — 1 
l^k — „ 


(19) 


By the additive property of variances under convolution, the discrete variances 
of the discrete temporal scale-space kernels will perfectly match those of the 
continuous model, whereas the mean values and the temporal delays may differ 
somewhat. If the temporal scale is large relative to the temporal sampling 
density, the discrete model should be a good approximation in this respect. 

By the time-recursive formulation of this temporal scale-space concept, the 
computations can be performed based on a compact temporal buffer over time, 
which contains the temporal scale-space representations at temporal scales 
and with no need for storing any additional temporal buffer of what has occurred 
in the past to perform the corresponding temporal operations. 


5.2 Discrete implementation of spatial Gaussian smoothing 

To implement the spatial Gaussian operation on discrete sampled data, we do 
hrst transform a spatial scale parameter in units of e.g. degrees of visual angle 
to a spatial variance s relative to a unit sampling density according to 

s=v^al ( 20 ) 

where p is the number of pixels per spatial unit e.g. in terms of degrees of visual 
angle at the image center. Then, we convolve the image data with the separable 
two-dimensional discrete analogue of the Gaussian kernel (Lindeberg |13j l 

r(ni,n 2 ; s) = (s)/^^(s) 


( 21 ) 
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where In denotes the modified Bessel functions of integer order and which cor¬ 
responds to the solution of the semi-discrete diffusion equation 


dsL{ni,n 2 ; s) = -{V\L){ni,n 2 ] s) 


( 22 ) 


where V^ denotes the five-point discrete Laplacian operator defined by (V^ /)(ni, n 2 ) 
f{ni - l,n 2 ) + /(ni -I- l,n 2 ) -I- f{ni,n 2 - 1) + /(ni,n 2 + 1) - 4/(ni,n2). These 
kernels constitute the natural way to define a scale-space concept for discrete 
signals corresponding to the Gaussian scale-space over a symmetric domain in 
the sense of guaranteeing non-enhancement of local extrema, while also ensuring 
a semi-group property T[-, •; si) * T{-, •; S 2 ) = T(-, •; si -I- S 2 ) over the discrete 
domain which implies that representations at coarser scales can be computed 
from representations at finer scales using a cascade property. 

This operation can be implemented either by explicit spatial convolution with 
spatially truncated kernels J2ni=-N 'l 2 n 2 =-N T{ni,n 2 ] s) > 1 — £ for small £ of 
the order 10“® to 10“® with mirroring at the image boundaries (adiabatic bound¬ 
ary conditions) or using the closed-form expression of the Fourier transform 

¥>(^ 1 ,^ 2 ) =E“=-ooE“=-oon^T^2; 


5.3 Discrete implementation of spatio-temporal receptive fields 

For separable spatio-temporal receptive fields, we implement the spatio-temporal 
smoothing operation by separable combination of the spatial and temporal scale- 
space concepts in sections |5.1| and |5.2| From this representation, spatio-temporal 
derivative approximations are then computed from difference operators 


— (—1, +1) 

d„ = (l,-2,l) 

(23) 

4 = (-2,0,+^) 

5,, = (1,-2,1) 

(24) 


expressed over the appropriate dimension. From the general theory in (Linde¬ 
berg nni) it follows that the scale-space properties for the original zero-order 
signal will be transferred to such derivative approximations, thereby implying 
theoretically well-founded implementation of discrete derivative approximations. 

For non-separable spatio-temporal receptive helds corresponding to a non¬ 
zero image velocity v = (ni,U 2 )^, we implement the spatio-temporal smoothing 
operation by first warping the video data (a;^,^^)^ = (a:i — vit,X 2 — V 2 t)'^ us¬ 
ing spline interpolation. Then, we apply separable spatio-temporal smoothing 
in the transformed domain and unwarp the result back to the original domain. 
Over a continuous domain, such an operation is equivalent to convolution with 
corresponding velocity-adapted spatio-temporal receptive helds, while being sig- 
nihcantly faster in a discrete implementation than corresponding explicit convo¬ 
lution with non-separable receptive helds over three dimensions. 

In addition to a transfer of the scale-space properties from the continuous 
model to the discrete implementation, all the components in this discretiza¬ 
tion, the discrete Gaussian kernel, the time-recursive hlters and the discrete 


Separable time-causal and time-recursive spatio-temporal receptive fields 


11 


derivative approximations, can be seen as mathematical approximations of the 
corresponding continuous counterparts. Thereby, the behaviour of the discrete 
implementation will approach the corresponding continuous model. 


6 Summary and discussion 


We have presented an improved computational model for spatio-temporal recep¬ 
tive fields based on time-causal and time-recursive spatio-temporal scale-space 
representation defined from a set of first-order integrators or truncated expo¬ 
nential filters coupled in cascade over the temporal domain in combination with 
a Gaussian scale-space concept over the spatial domain. This model can be ef¬ 
ficiently implemented in terms of recursive filters and we have shown how the 
continuous model can be transferred to a discrete implementation while retain¬ 
ing discrete scale-space properties. Specifically, we have analysed how remain¬ 
ing design parameters within the theory, in terms of the number of first-order 
integrators coupled in cascade and a distribution parameter of a logarithmic 
distribution, affect the temporal response dynamics in terms of temporal delays. 

Compared to other spatial and temporal scale-space representations based 
on continuous scale parameters, a conceptual difference with the temporal scale- 
space representation underlying the proposed spatio-temporal receptive fields, 
is that the temporal scale levels have to be discrete. Thereby, we sacrifice scale 
invariance as resulting from Gaussian scale-space concepts based on causality or 
non-enhancement of local extrema (Koenderink [lb] ; Lindeberg m) or used as 
a scale-space axiom in certain axiomatic scale-space formulations (lijima |17| : 
Florack et al. [15]; Pauwels et al. [151; Weickert et al. [55]; Duits et al. [5T] : 
Fagerstrom [5^); see also Koenderink and van Doom [53], Florack et al. [53] 
and Tschirsich and Kuijper [^ for other scale-space approaches closely related 
to this work. For a vision system intended to operate in real time using no other 
explicit storage of visual data from the past than a compact time-recursive buffer 
of spatio-temporal scale-space at different temporal scales, the loss of a contin¬ 
uous temporal scale parameter may however be less of a practical constraint, 
since one would anyway have to discretize the temporal scale levels in advance 
to be able to register the image data to perform any computations at all. 

In the special case when all the time constants of the first-order integrators 


are equal, the resulting temporal smoothing kernels in the continuous model (141 


correspond to Laguerre functions, which have been previously used for modelling 
the temporal response properties of neurons in the visual system (den Brinker 
and Roufs [55]) and for computing spatio-temporal image features in computer 
vision (Rivero-Moreno and Bres [57]; Berg et al. [55]). Regarding the corre¬ 
sponding discrete model with all time constants equal, the corresponding dis¬ 
crete temporal smoothing kernels approach Poisson kernels when the number of 
temporal smoothing steps increases while keeping the variance of the composed 
kernel fixed (Lindeberg and Fagerstrom [35] )• Such Poisson kernels have also 
been used for modelling biological vision (Fourtes and Hodgkin [55]). Gompared 
to the special case with all time constants equal, a logarithmic distribution of 
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the intermediate temporal scale levels (|^ does on the other hand allow for larger 
flexibility in the trade-off between temporal smoothing and temporal response 
characteristics, specifically enabling faster temporal responses (shorter temporal 
delays) and higher computational efficiency when computing multiple temporal 
or spatio-temporal receptive field responses involving coarser temporal scales. 
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